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We consider the two-layer Heisenberg antiferromagnet near a zero temperature quantum phase 
transition from a disordered dimer phase to a collinear Neel state. At approaching the transition 
point the spin-wave gap vanishes as A a (Jj_ — J± c ) v . To account for strong correlations between 
the S = 1 elementary excitations we apply the Brueckner diagram approach which gives the critical 
index v ~ 0.5. We demonstrate also that the linearized in density Brueckner equations give the 
mean field result v = 1. Finally an expansion of the Brueckner equations in powers of the density, 
combined with the scaling hypothesis, give v ~ 0.67. This value reasonably agrees with that of 
the nonlinear 0(3) a-model. Our approach demonstrates that for other quantum spin models the 
critical index can be different from that in the nonlinear cr-model. We discuss the conditions for this 
to occur. 
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There is a significant interest in the physics of quan- 
tum phase transitions in S = 1/2 two-dimensional 
(2D) Heisenberg models in connection with the high-T c 
cuprate superconductors. It has been argued by Hal- 
dane |l| and Chakravarty, Halpcrin and Nelson Q that 
the phase transition is described by the 2+1-dimcnsional 
non-linear 0(3) a model predicting that at approach- 
ing the transition point the spin-wave gap vanishes as 
A ~ (J_l — J± c ) 1 ' with critical index v sa 0.7. Assum- 
ing the cr-model description the universal dynamic and 
static properties of two-dimensional antiferromagnets in 
the vicinity of a zero-temperature phase transition have 
been studied in detail by Sachdev, Ye and Chubukov 
Hfj}. However mapping of the Heisenberg model to the 
cr-model has been proven only for the one-layer square 
lattice antiferromagnet which itself does not exhibit a 
zero temperature quantum phase transition. 

In the present work we investigate the two-layer 
Heisenberg antiferromagnet which has a zero tempera- 
ture quantum phase transition from a disordered dimer 
phase to a collinear Neel phase as the dimerization de- 
creases. For this model the numerical results of Singh, 
Gelfand and Huse || and Weihong || obtained by se- 
ries expansions, as well as Monte Carlo calculations per- 
formed by Sandvick and Scalapino |?j demonstrate rea- 
sonable consistency of the critical index with that in the 
cr-model. At the same time Sandvik and Vekic || us- 
ing quantum Monte-Carlo simulations have shown that 
for some other 2D antiferromagnetic models (dimers ar- 
ranged in ladders and dimers in a staggered pattern) 
the critical exponents are different from the nonlinear 
cr-model prediction. In this situation it is very important 
to analyze the behavior near the critical points by an in- 
dependent analytical method. Such a method, based on 
the Brueckner perturbation theory, has been developed 
in the paper M . In the present work we apply it for the 
calculation of the spin- wave gap critical index in the two- 
layer Heisenberg antiferromagnet. This method gives an 
independent calculation of the index for this model, and 
additionally it allows us to formulate the conditions when 
one could expect deviations from the cr-model behavior. 

We consider the two-layer, S = 1/2, Heisenberg anti- 
ferromagnet described by 

= J ]T (SiiS w + S 2i S 2j ) + Jl s«s 2i . (1) 

<i,j> i 

The spins Si and S2 represent two planes coupled by 
J_i_. Both couplings are antiferromagnetic (J,J± > 0) 
and the first sum runs over nearest neighbors on a square 
lattice. Considering the "J-terms" in (|l|) as a perturba- 
tion one can introduce the bond operator representation 
^0| , p| and exactly map the Hamiltonian (|l|) to the effec- 
tive Hamiltonian ||: 
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where t t is the creation operator of the triplet at 
the bond i and a = x, y, z is the polarization of the 
triplet. The operator is the Fourier transform of 

*L : *k,a = ^E r *r, Q e l ( k+k °) r . As usual the momen- 
tum takes values inside the Brillouin zone — 7r < k x < tt, 
— 7r < k x < 7T, but we shift the argument in the Fourier 
transform by ko = (tt, it). In this notation the mini- 
mum of the spin- wave dispersion is at k = 0. The co- 
efficients are of the form = J± + 2J£k, B^ = 2 J£k, 
£k = — ^{cosk x + cosfcy). An infinite on-site repulsion 
between triplets Hjj is introduced to take into account 
the hard-core constraint = (only one triplet can 

be excited on a bond). The interaction Hy gives the 
dominant contribution to the renormalization of the spin- 
wave spectrum. It has been demonstrated |J that in 
the Brueckner approximation this renormalization is de- 
scribed by the self-energy operator 



S(k,w) 



-5> q <r(k + q, w 
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where the scattering amplitude is 




(7) 



and the spin wave spectrum Wk, quasiparticle residue Z k 
and Bogoliubov parameters Uk, fk are given by the for- 
mulas: 



R2 



(8) 



I k = J x + 2 J^k + E(k, 0) + 4J£ k J2 ti z i v l, (9) 

q 



H e ff — H2 + H4 + Hi 



(2) 
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These equations take also into account the quartic inter- 
action (pj) in one-loop approximation. In order to find 
spectrum, equations (U0|||9|;|lO|) have to be solved self- 
consistently for S(k, 0) and Z^. The plot of the spin- wave 
gap A versus J±/J is presented in Fig.l. In the same fig- 
ure we present results from the papers (6|^] obtained by 
the dimer series expansions. 

Close to the critical point (A <C J) and for small mo- 
menta (fc <C 1) the dispersion can be represented as 



uj k w y/A 2 + c 2 k 2 , 



(11) 



where c = 1.85 J is the spin- wave velocity §|||l§. To 
find the critical index let us write eq. m) at the point 
k = 



(12) 



It is convenient to introduce the values of Aq and Bq at 
the critical point: Aq c = —Bq c . Let us vary Jj_ keeping 
J fixed, and let us introduce the deviation from the crit- 
ical point 5J± = Jj_ — J_\_ c . Considering SJ± and A as 
independent variables and using eqs. (|9| JlO| , pT|) we find 
the variations of Aq and Bq : 



A = A 0c + 5J± + 5E(0, 0) - A, 



(13) 



B = B Qc + A. 



To find the variation of the self-energy notice that accord- 
ing to eq. (^) T(q, — w q ) ocgatA/c~g<Cl. Therefore 
from eq. (fj) we find 

J 



6E(0, 0) oc 1 5v \ oc A 2 In - 



(14) 



where x = J±/J and A is found from eqs. (|6j]7|,|8|,|9|,[l0|). 

For pure scaling behavior [A oc (8J^) V ] the function 
fi(x) is linear and the function ji gives the index: 
f2(x) = v. Keeping in mind that when the gap is small 
(ln(J/A) 3> 1) our approximation is not justified, we con- 
clude from the plot of /2 that the estimate of the critical 
index is v G [0.6,0.75] or 



0.67 ±0.07. 



(17) 



Another way to look at the critical index is to expand 
the Brueckner equations in powers of the density rtj,. In 
the leading approximation Z q = u q = 1 and therefore 
the vertex (|t]) takes the form 



r(q, -w q ) 



d 2 p 



(27r) 2 W q + UJ p + ( 



"p-q 



(18) 



In contrast to ([?]) it does not vanish at q — » 0. Let us de- 
note r(0, 0) = F c . The last terms in the equations jl^ ) 
are due to the quartic interaction H4 (|J). This inter- 
action is relatively small and therefore we neglect these 
terms. Below we prove that they are really small. So 
instead of eqs. (|l^) we have 



A = A 0c + SJ X + 5E(0, 0) 

Bq = B 0c . 



(19) 



If we substitute this into eq.([L2|) and neglect terms 
quadratic in A we find that variation of Aq must van- 
ish 

SA = 8J ± + 6E(0, 0) = 0. (20) 
The variation of the self-energy should be found from eq. 



Terms linear in A are canceled out after substitution of 
( |l3| ) and ([l4|) into eq. (|l2]), and neglecting logarithmic 
dependence we find 



A oc 



(15) 



Thus the critical index in the Brueckner approximation is 
v = 0.5. It is known from nuclear and atomic physics |l5| ] 
that the Brueckner approximation usually works pretty 
well even in systems with high density, but parametri- 
cally it is justified only at the low density limit. In essence 
it is the dilute gas approximation. It has been demon- 
strated m that for the model under consideration the 
actual small parameter is nf,ln(J/A), where rib ~ 0.1 is 
density of the triplet excitations. It is clear that when the 
gap is very small this parameter is becoming large and 
the gas approximation can fail. To analyze the situation 
numerically we plot in the Fig. 2 the functions 



AW 



dx ) 



(16) 



, . , ( ^ln(A/J) \ 2 /rf 2 ln(A/J) 



5S(0,0)=4| ^r(q,- Wq )^ 2 +4 1 ^«5r(q,-u; q )<. 

(21) 

The main contribution to the first integral in this formula 
comes from small momenta (q ~ A/c <C 1) since 



0J a \ LO a 




A 0c A 2 



4(A 2 +cV) 3/2 



d 2 s 



(22) 

The variation of A in this formula vanishes according to 
eq. (^0|). Then the integral can be easily evaluated and 
the result is — r c Ao c A/7rc 2 . The main contribution to 
the second integral in ( plf ) comes from large momenta 
(q ~ 1) where we can write 8T(q, —ojq) = T'8J±. It is 
obvious from ([l8]) that at J± ^ J the derivative V = 3. 
However we need this derivative near the critical point 
where the numerical calculation shows that T' 2.9 in- 
dependent of momenta. Altogether the variation ( pl| ) can 
be represented as 



3 



££(0,0) 



r c A 0c 



A- 



-T'ni,6J± 
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where rib = Sj^^Vq ~ 0-12 is the density of spin- wave 
excitations at the critical point. At this step we can check 
how small is the neglected quartic interaction term in eqs. 
df). The ratio of this term to the "A-term" in ( pi) i s 
J/r c « 0.15, since T c w 6.3 J. After substitution of (|2j) 
into eq. ( ppf ) we find the relation between the gap and 
5J X 



TTC 4 

——(1 + -r'n 6 )(J_L - Jxc) « - Ji 

1 c^Oc o 



(24) 



The coefficient 1.1 corresponds to r c w 6.3J, Aq c rs 2.4 J, 
c w 1.9J, r' w 2.9, and rib w 0.12 which have been found 
by numerical solution of the linearized in density eqs. 
(p|,pPJlO|). We remind that the linearization means that 
the vertex is taken from eq. ( |l8| ) and the residues in (|^) 
are replaced by unity (Z q = 1). The numerical solution 
also gives A(J±) shown in Fig.l. The slope at J± = J± c 
is in perfect agreement with the semianalytical formula 
(p4|). Thus the leading term in powers of the density 
gives the critical exponent v = 1, the same as the mean 
field approximation. 

Now, consider the first correction due to the triplet 
density rib. We will keep only those terms which con- 
tain additional lnc/ or ln(A/J), so the parameter is 
rib ln(A/ J). The rib ln(A/J) terms arise only from expan- 



sion of the vertex (Q), replacing 



t>q_ p wc obtain 



p q-p 



T(q, 



d 2 p 1 



(2tt) 2 LU q 



- w c 



= (l+vi)(l + 



(25) 
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After simple integration one can find that for small q 
(?~ A/c«l) the vertex is 



r(g,-u; g ) 



Anc 2 



liiq. 



(26) 



Substitution into (21) gives the variation of the self- 
energy 



< 5E(0,0)«-^ £ A| 1 



r A, A , 4 . 



7TC 



In— + -rn b SJ ± , 



4ttc 2 J / 3 



(27) 



which together with eq. ( |20| ) results in 

4„, ... A r c i 0c , 5J± 



A 



7TC 



- -(1 + -T'n b )SJ± 1 



In- 



Let us assume the scaling behavior A oc (SJ±) 1 ' with 
v = 1 — /3. Expanding this formula in powers of /3 and 
comparing with (Eq) we find 



v = 1 - 



47TC 2 



0.67, 



(29) 



which agrees with the estimate (|17j) a nd with the result 
of the cr-model approach(z/ sa 0.70 |16|). 

In conclusion, using the Bruekner approach for two- 
layer Heisenberg antiferromagnet we have calculated the 
spin-wave critical index for the zero temperature quan- 
tum phase transition. The result is in reasonable agree- 
ment with that of the nonlinear 0(3) cr-model. This 
agreement is due to the relative smallness of the quar- 
tic interaction for the model under consideration. In this 
situation the hard-core constraint is the most important 
and it is very natural that the result is similar to that 
of the cr-model. However this situation is not general. 
There are many models where the quartic interaction is 
very important. It can even produce bound states of 
triplet spin waves |14j which effectively change the num- 
ber of relevant degrees of freedom. In this situation one 
can expect a very substantial deviation from the simple 
cr-model. An important example of such a system is the 
2D Ji — J2 model where the singlet bound state has an 
extremely low energy [[[5) . 

We are very grateful to V. N. Kotov, R. R. P. Singh 
and Z. Weihong for stimulating discussions. This work 
was supported by a grant from the Australian Research 
Council. 



[1 
[2 

[3 
[4. 

[5 

[6 

ir 

[8] 

[9] 

[10 

[11 
[12 



(28) [13] 



e-mail: pavel@newt.phys.unsw.edu.au 
F.D.M.Haldane, Phys.Rev.Lett. 50, 1153 (1983) 
S.Chakravarty, B.I.Halperin, and D.R.Nelson, 
Phys.Rev.Lett. 60, 1057 (1988) 
S.Sachdev, J.Ye, Phys.Rev.Lett. 69, 2411 (1992) 
A.V.Chubukov, S.Sachdev, and J.Ye, Phys.Rev.B 49, 
11919 (1994) 

R.R.P.Singh, M.P.Gelfand, and D.A.Huse, 
Phys.Rev.Lett. 61, 2484 (1988) 
Z.Weihong, Phys.Rev.B 55, 12267 (1997) 
A.W.Sandvik and D.J.Scalapino, Phys.Rev.Lett. 72, 
2777 (1994) 

A.W.Sandvik and M.Vekic, J. of Low Temp.Phys. 99, 
314 (1995) 

V.N. Kotov, O.Sushkov, Z.Weihong and J.Oitmaa, Phys. 

Rev. Lett. 80, 5790 (1998). 

A. V. Chubukov, JETP Lett. 49, 129 (1989). 

S. Sachdev and R. Bhatt, Phys. Rev. B 41, 9323 (1990). 

K.Hida, J.Phys.Soc.Jpn 61, 1013 (1992); M.P.Gelfand, 

Phys.Rev.B 53, 11309 (1996) 

V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov and O. 
P. Sushkov. "Relativistic Quantum Electrodynamics and 
Weak Interaction Effects in Atoms" AIP conference pro- 
ceedings, 189, p 254 (1989). Editors: W. Johnson, P. 
Mohr, J. Sucher. 



4 



[14] O. P. Sushkov and V. N. Kotov. Phys. Rev. Lett. 81, 
1941 (1998). 

[15] V. N. Kotov, J. Oitmaa, O. Sushkov, and Z. Weihong, 

work in progress. 
[16] M.Ferer and A.Hamid-Aidinejad, Phys.Rev.B 34, 6481 

(1986), and references cited therein 

FIGURE CAPTIONS 

Fig.l Triplet gap A/ J as a function of x = J±/J- Solid 
and dashed lines are the results of the self-consistent so- 
lution using the Bruekner approach and the Bruckner 
equations linearized in density respectively. The dots 
(with error bars) are estimates obtained by dimer series 
expansions. 

and f 2 (x) = 
/ d lD }^^ J ^ versus x = Jj_/J (sec text). 



Fig. 2 Functions f\(x) 
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